Super-diffusion around the rigidity transition: 
Levy and the Lilliputians 
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By analyzing the displacement statistics of an assembly of horizontally vibrated bidisperse fric- 
tional grains in the vicinity of the jamming transition experimentally studied before [l|, we establish 
that their superdiffusive motion is a genuine Levy flight, but with 'jump' size very small compared 
to the diameter of the grains. The vibration induces a broad distribution of jumps that are random 
in time, but correlated in space, and that can be interpreted as micro-crack events at all scales. As 
the volume fraction departs from the critical jamming density, this distribution is truncated at a 
smaller and smaller jump size, inducing a crossover towards standard diffusive motion at long times. 
This interpretation contrasts with the idea of temporally persistent, spatially correlated currents 
and raises new issues regarding the analysis of the dynamics in terms of vibrational modes. 

PACS numbers: 



INTRODUCTION 



As the volume fraction of hard grains is increased be- 
yond a certain point, the system jams and is able to sus- 
tain mechanical stresses. This rigidity /jamming transi- 
tion has recently been experimentally investigated [l], Q 
in an assembly of horizontally vibrated bi-disperse hard 
disks, using a quench protocol that produces very dense 
configurations, with packing fractions beyond the glass 
density <p g , such that the structural relaxation time r a is 
much larger than the experimental time scale. There is 
a density range <f> Q < <p < qb a where the strong vibration 
can still induce micro-rearrangements through collective 
contact slips that lead to partial exploration of the por- 
tion of phase space restricted to a particular frozen struc- 
ture. For (j) g < (j) < 4>j ~ 0.842, the system is frozen but 
not rigid; the system can only sustain an external stress 
for cj) larger than 0j, which appears as a genuine critical 
point where a dynamical correlation length and a cor- 
relation time simultaneously diverge, showing that the 
dynamics occurs by involving progressively more collec- 
tive rearrangements. 

One of the most surprising result of [l[ was the dis- 
covery of a superdiffusive regime in the vicinity of qbj. 
In a time range that diverges when <f> — >> 0j, the typical 
displacement of the grains grows as r u , with v > 1/2, 



i.e. faster than the familiar diffusive yfr law, while al- 
ways remaining small compared to the diameter of the 
grains! It was also found that dynamical heterogeneities 
reach a maximum when r = t* ((/>), precisely when the 
superdiffusive character of the motion is strongest Q , in 
agreement with the general bound on the dynamical sus- 
ceptibility established in [3, 4]. This superdiffusion was 
interpreted as the existence of large scale convective cur- 
rents, which were tentatively associated to the extended 
soft modes that appear when the system loses or acquires 
rigidity at |7], and along which the system should 
fail. A similar interpretation of the dynamics of partic- 
ulate systems close to the glass or jamming transitions 
was promoted in [8]. The intuitive idea is that energy 
barriers that the system has to cross to equilibrate are 
essentially in the directions (in phase space) defined by 
the soft modes. Within this picture, the harmonic motion 
of the particles around a metastable configurations can 
be used to guess the structure of the anharmonic barrier 
crossing events. This is certainly reasonable when the 
barriers are small, so that the top of the barrier is still in 
a quasi-harmonic regime. 

The primary aim of the present paper is to revisit the 
above interpretation for our system of hard-disks with 
friction, in the light of a deeper analysis of the statistics 
of displacements. To our surprise, we find that the distri- 
bution of rescaled displacements is, in the superdiffusive 
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FIG. 1: Left: Sketch of the experimental setup. Right: Pic- 
ture of the grains together with an illustration of the relax- 
ation pattern of the field qi at r* (see text below for details). 



regime, accurately given by an isotropic Levy stable dis- 
tribution L M . Furthermore, the exponent fi characteriz- 
ing this Levy stable distribution is equal to 1/^, where v 
is the superdifTusion exponent, precisely what one would 
expect if the motion of the grains was a Levy flight, 
i.e. a sum of uncorrelated individual displacements with 
a power-law tail distribution of sizes such that the vari- 
ance of the distribution diverges [H, ^ . This divergence 
only occurs at <p = but is truncated away from the 
critical point, which explains why the motion reverts to 
normal diffusion at very large times. This finding shows 
that the rearrangements corresponding to the maximum 
of dynamical correlations cannot be thought of as large 
scale currents that remain coherent over a long time scale 
r* . SuperdifTusion is not induced by long-range temporal 
correlations of the velocity field, as was surmised in |l|. 
Quite on the contrary, the total displacement on scale 
r* is made of a large number of temporally incoherent 
jumps with a broad distribution of jump sizes. As is well 
known, a Levy flight is, over any time r, dominated by a 
handful of particularly large events. Hence, the only pre- 
dictability of the total displacement between t = and 
t = r* on the basis of the motion in any early period, 
say between t = — >• r*/10, comes from the presence of 
one of these large jumps in [0,r*/10]. The situation for 
frictional grains might therefore be quite different from 
what happens in supercooled liquids or frictionless hard 
spheres, when soft modes should indeed play an impor- 
tant role. 

On the other hand, the very fact that individual jumps 
are broadly distributed in such a dense system where 
grains can hardly move means that these large jumps are 
necessarily collective, with a direct relation between the 
size of the jump and the number of particles involved that 
we discuss below. The physical connection between su- 
perdifTusion and dynamical heterogeneities in this system 
becomes quite transparent. Note however that although 
we speak about "broad distributions" and "large jumps", 
one should bear in mind that all of this takes place on 
minuscule displacement scales, a few 10 ~ 2 of the grains 
diameter! Our Levy flights are therefore Lilliputian walks 
with a diverging second moment... 

The experimental set-up and the quench protocols are 
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FIG. 2: Left: Typical trajectories of a grain during a 10000 
cycle acquisition. The packing fractions are (f) — 0.8402 
(red), 0.8417 (orange), 0.8422 (yellow), 0.8430 (green). 
would stand between the orange and yellow trajectories. 
Right: Fraction of broken neighborhood relations (links in 
the Voronoi' tessellation sense), as a function of time for all 
packing fractions. 

described in detail in [1], and summarized in Fig. [TJ 
The stroboscopic motion of a set of 8500 brass cylin- 
ders ("grains") is recorded by a digital video camera. 
The cylinders have diameters d\ — 4 ± 0.01mm and 
d,2 = 5 ± 0.01mm and are laid out on a horizontal glass 
plate that harmonically oscillates in one direction at a 
frequency of 10 Hz and with a peak-to-peak amplitude 
of 10 mm. The cell has width L w 100 d sma u, and its 
length can be adjusted by a lateral mobile wall controlled 
by a fim accuracy translation stage, which allows us to 
vary the packing fraction of the grains assembly by tiny 
amounts ~ 5 10 -4 ). The position of the grains is 

tracked within an accuracy of2.10~ 3 di. In the following, 
lengths are measured in d\ units and time in cycle units. 

STATISTICAL CHARACTERIZATION OF THE 
MOTION OF GRAINS 

Fig. [2]-left shows the trajectories of a single grain dur- 
ing the whole experimental window for packing fractions 
sitting on both sides of the transition <pj. As expected, 
the typical displacements are much larger at lower pack- 
ing fractions. The so-called caging dynamics appears 
clearly on the trajectories obtained at larger packing frac- 
tions: the grains seem to stay confined around fixed posi- 
tions for a long time and then hop around from one cage 
to another over longer time scales. The striking feature 
of these curves is the remarkably small amplitude of the 
motion that we alluded to above. Even at the lowest 
packing fractions, the grains do not move much further 
than a fraction of a diameter over the total course of one 
experimental run. Another evidence of the absence of 
structural relaxation is given by the low percentage of 
neighboring links (in the Voronoi sense) that are broken 
during the 10 4 cycles of an experimental run (see Fig. [2} 
right). This fraction goes from around 5% for the loosest 
packing fraction to around 0.2% for the densest. This 
confirms that for densities around </>j, the system is deep 
in the glass phase where structural relaxation is absent, 
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FIG. 3: Top: Average displacement (r) = (|r-|)^t (left); 
and local slope v — dlna^ (r)/d\nr (right) as functions of 
the lag r for all studied values of cj). The dotted dark line 
corresponds to <j>j. Bottom left: Characteristic time scales 
extracted from the analysis of v(t) including r*, which will be 
introduced and discussed in the section devoted to dynamical 
heterogeneities. Bottom left: Pdf of the displacements n, 
P r (r), normalized by 0>(r), for different values of the lag 
time r at cj) = 0.8446 (right). A Gaussian distribution (black 
dotted curve) is shown for comparison. 

i.e. in a very different regime from the one studied in 

§1111. 

We will denote by = (x^yj) the vector position of 
grain i at time t in the center-of-mass frame. Lagged dis- 
placements are denned as y\{t) = R^ +r — R^ for grain i 
between time t and t-\-r. We have checked in details that 
the statistics of v\ along the x and the y axes are identi- 
cal, or more precisely that the motion is isotropic, in spite 
of the strongly anisotropic nature of the external drive. 
This in itself is a non trivial observation, that shows that 
the random structure of the packing is enough to convert 
a directional large scale forcing into an isotropic noise 
on small scales. From now on, we will thus focus on 
the total displacement r\ = |r-|. We measure typical dis- 
placements for a density (j) and lag r as the mean absolute 
displacement, denned as: 

MT) = <|r|(r)|>. t (1) 

where the time average (-) t is performed over 10,000 cy- 
cles. We have chosen to estimate cr^r) by the mean 
absolute value deviation instead of the root mean square 
displacement because - as we shall see - the later diverges 
when approaching the transition. 

The evolution of cr^(r) with the packing fraction is pre- 



sented in Fig. [3j see Note the very small values of 
cT(f)(r) <C di at all timescales, which is a further indica- 
tion that the packing indeed remains in a given struc- 
tural arrangement. At low packing fractions (j) < <pj and 
at small r, cr^(r) displays a sub-diffusive behavior. At 
longer time, diffusive motion is recovered. As the packing 
fraction is increased, the typical lag at which this cross- 
over occurs becomes larger and, at first sight, cr^(r) does 
not seem to exhibit any special feature for <\> = </>j (corre- 
sponding to the bold line in Figj3j). Above </>j, an inter- 
mediate plateau appears before diffusion resumes. How- 
ever, a closer inspection of cr^r) reveals an intriguing 
behavior, that appears more clearly on the local logarith- 
mic slope v = <9 In <T0(r)/d In (t) shown in Fig. [3j When 
v = ^, the motion is diffusive, whereas at small times, 
v < \ indicates sub-diffusive behavior. At intermedi- 
ate packing fractions, instead of reaching \ from below, 
v overshoots and reaches values ^0.65 before reverting 
to ^ from above at long times. Physically, this means 
that after the sub-diffusive regime commonly observed 
in glassy systems, the particles become super- diffusive 
at intermediate times before eventually entering the long 
time diffusive regime. At higher packing fractions, this 
unusual intermediate super diffusion disappears: one only 
observes the standard crossover between a plateau regime 
at early times and diffusion at long times. In order 
to characterize these different regimes, we define three 
characteristic times: ti(0) as the lag at which v{r) first 
reaches 1/2, corresponding to the beginning of the super- 
diffusive regime, r s jj((j)) when v(t) reaches a maximum 
z/ *(0) (peak of super-diffusive regime), and rjj((j)) where 
v{t) has an inflection point, beyond which the system ap- 
proaches the diffusive regime. These characteristic time 
scales are plotted as a function of the packing fraction in 
the bottom-left panel of Fig. [3l Whereas t\ does not ex- 
hibit any special feature across (/>j, both t s d and td are 
strongly peaked at <pj. We have also shown on the same 
graph the time scale r* where dynamical heterogeneities 
are strongest (see [1] and below). As noted in the intro- 
duction, r* and r s jj are very close to one another and we 
will identify these two time scales in the following. 

We now turn to the probability distribution of the dis- 
placements r\ (r). We have represented on the bottom- 
right panel of Fig. [3] the distributions accumulated for 
all grains and instants for the packing fractions closest 
to and several values of lag-time r. The horizontal 
axis has been normalized by the root mean square dis- 
placements, and a unit Gaussian is also plotted for com- 
parison. We find that the tails of the distributions are 
much fatter than the Gaussian, unveiling the existence of 
extremely large displacements (compared to typical val- 
ues) in the region of time lags corresponding to superdif- 
fusion. As we will discuss below, these fat tails tend 
to disappear when one leaves the super diffusion regime, 
i.e. when \<j> — <j) c \ and/or | lnr/r*| become large. 

Before characterizing more precisely these probability 
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distributions, we would like to mention that we have re- 
moved "rattling events" from our analysis, i.e. events 
where particles make an anomalously large back and 
forth motion, of amplitude larger than (say) O.ldi. Some 
of these events are due to the same grain rattling during a 
large fraction of the experimental run, while other events 
are localized in time for a given grain, and could actually 
well be part of the extreme tail of the distribution seen 
in Fig. 02 Our point here is that the observed fat-tails 
are not an artifact due to a few "loose" grains - each and 
everyone of the grains seem to contribute at one point or 
another in its history to these fat tails. 

A neat way to characterize the probability distribution 
of the displacements is to study the generating function 
of the squared displacements: 



k=0 



<rP(r)). (2) 



As a benchmark, one can easily compute T (A, r) in the 
case of isotropic Gaussian diffusion. One obtains: 



1 



1 



(3) 



where x = 2(j0(r) 2 A is the scaled variable. When com- 
puting T (A,r) for the empirical data, one finds, as ex- 
pected, systematic discrepancies with the T (A,r) G . The 
rescaling of the different curves when plotted as a func- 
tion of x, on the other hand, works very well, as can be 
inferred from the scaling of the probability distributions 
themselves. 

The most important finding is that the small x be- 
havior of the empirical ^(X^r) is singular: as shown in 
the top-right panel of Fig. 3J l/J r (A,r) — 1 behaves as 
x a when x <C 1, with a ~ 0.80. It is easy to show 
that such a singular behavior is tantamount to the ex- 
istence of a power-law tail in the distribution of r(r), 
P T (r) ~ r ^oc r _1_M with a = min(l,/i/2). The value of 
a therefore corresponds to fi « 1.6 < 2, which means 
that the variance of the distribution of r is formally di- 
vergent, or at least dominated by a large physical cut- 
off r max . This suggests to fit ^(X^r) over the whole x 
regime by the appropriate Laplace transform of a Levy 
stable distribution L M , that reads for the present isotropic 
two-dimensional problem [12| : 



F{\,t) l =F l {x) 



f 

Jo 



dz e 



-C(xz) 



A'/ 2 



(4) 



where fi < 2 and C is a numerical constant. It is easy to 
check that for small x, T L (x) « 1 - CT(1 + f ) x^ 2 . Our 
main result is that this function is a very good fit to the 
empirical data corresponding to <p « <pj and r « r*, for 
all values of x, see the bottom-left panel in Fig. HJ The 
optimal value of fi is slightly smaller, ji « 1.6, than the 
one obtained by a fit of the small x region. The quality of 
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FIG. 4: Top: T T {x) as a function of — log 10 (x) for r = 
r* = 373 (dots, left panel) and l/J 7 T (x) — 1 as a function of 
x for different values of r (dots, right panel), from 1 (blue) 
to 10 4 (red). For both plots: (j) = (f>j and the solid line 
corresponds to the Gaussian fit 1/(1 + Bottom-left: Fit of 
T T (x) for (f) & cj)j and r ~ r* using a Levy distribution model 
J 7 (A,t) l with index ji = 1.6, and comparison with a Student 
distribution model T(X,t) l with the same small x behavior. 
Bottom-right: Fitting parameter s = 1 — /i/2 as a function of 
6 and r. 



the fit suggests that the distribution of displacements is 
indeed an isotropic Levy stable distribution L M of index 
fi. 

We have checked that other distribution functions with 
power-law tails achieve a much poorer fit to the data. For 
example, a Student distribution of displacements (which 
cannot be obtained as the sum of independent jumps, 
contrarily to the Levy distribution) leads to: 



= i - (c'xy' 2 



(5) 



(z + Cxy i 2 ' 



where \i < 2 and C is another numerical constant. Al- 
though T(\,t) s has the same small x singularity and 
the same large x behavior as J r (A,r) L , it turns out not 
to be possible to adjust the parameter C in such a way 
to adjust simultaneously the small and large x behavior 
of the empirical T{\,t). 

What is most significant, however, is that the value of 
the Levy index \i corresponding to the best fit turns out 
to be very close to the one expected in the case of a Levy 
flight where the diffusion exponent v is given by l//i, 
since we find v w 0.65 (see Fig. [3]) whereas l//i « 0.625. 

Fig. H] bottom-right gives the value of the fitting pa- 
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rameter e — 1 — /jl/2 in the plane (/>, r, obtained from the 
small x behavior of ^(A, r), corresponding to large dis- 
placements. We see that the effective values of e (resp. 
/i) become closer to (resp. 2), corresponding to normal 
diffusion, as when \<j> — <j) c \ and/or |lnr/r*| increase. We 
believe that this corresponds to a so-called "truncated" 
Levy flight [HI , with a cut-off value in the distribution of 
elementary jump size that becomes smaller and smaller 
as one departs from the critical point, rather than a con- 
tinuously varying exponent /i. In other words, a plausible 
scenario that explains the behavior of e as a function of 
0, r is that the tail of the elementary jump size distribu- 
tion is given by: 



Pi{r) 



(6) 



where the exponent /i ~ 1.6 is independent of (/>, whereas 
the typical scale of the jumps r*o(</>) and the cut-off 
^max(^) both depend on <\>. It is reasonable that r max ((/>) 
only diverges at <fi = 0j, while ro(</>) has a regular, de- 
creasing behavior as a function of <p. In fact, the parame- 
ter C appearing in the Levy distribution above is directly 
proportional to r£. 

Assuming that the jumps are independent, the distri- 
bution of displacements P T (r) on a time scale r is given 
by the r-th convolution of -Pi(r), which converges to a 
Levy distribution of order \i when r max = oo and pre- 



dicts r (xr v with v — l//i (see e.g. [6]). For finite r max , 
P T will be very close to in the intermediate regime 
1 « t « td, where r D ~ (r^^/r^Y, before crossing 
over to a diffusive regime at very long times. If r max (</>) 
behaves as \<p — </>j| - ^, the diffusion time td should di- 
verge as \<j) — </>j|~^ M when approaching <\>j. 

In order to check directly whether the jumps are in- 
deed independent, we have measured the correlation of 
instantaneous (r = 1) velocities (r-.r- +t )^. As shown 
in Fig. [5j this correlation function indeed decays ex- 
tremely fast with the lag t' \ excluding that superdiffusion 
could be due to long-range correlations in the displace- 
ments. However, we discovered that the amplitude of 
the displacements, proportional to the parameter 7*0, re- 
veal long-range correlations, decaying approximately as 
— lnt f . This suggests that some slow evolution takes 
place during an experimental run, that affects the value 
of ro over time. A natural mechanism is through the 
fluctuations of the local density <p (see below), that feed- 
back on 7*0. The approximately logarithmic decay of the 
correlations is interesting in itself and characteristic of 
multi-time scale glassy relaxation. These long-ranged 
correlations slow down but fortunately do not jeopardize 
the convergence towards a Levy stable distribution. [17] 
Therefore, the above interpretation in terms of a jump 
dominated superdiffusion, instead of persistent currents, 
is warranted. 




FIG. 5: Correlation of the instantaneous (r = 1) velocity 
vectors |(r*.r* +t )ij\ (blue circles) versus the correlation of 
their amplitudes (r*.r* +t )i,t (red squares) as functions of the 
lag t\ for cj) « (f)j. 



DYNAMICAL CORRELATIONS 

The above analysis focused on single grain statistics 
and established that the displacement of a grain has 
a Levy stable distribution in the superdiffusive regime. 
This means that the motion of each grain can be de- 
composed in a succession of power-law distributed jumps. 
But if "large" jumps can occur in such a jammed system, 
this necessarily implies that these jumps are correlated 
in space. In order to characterize these dynamical cor- 
relations, we have introduced in [1] the following local 
correlation function: 



„t2 



q\ (a, r) = exp - 



(r) 



2a 2 



(7) 



where a is a variable length scale over which we probe the 
motion. Essentially, qj (a, r) = if within the time lag r 
the particle has moved more than a, and qj (a, r) = 1 oth- 
erwise: q\ (a, r) is a mobility indicator. The average over 
all i and t defines a function akin to the self-intermediate 
scattering function in liquids: 



»(a,r) 



;^^^(a,r)) £ . 



(8) 



Note that it is closely related to the previously introduced 
generating function through Q(a,r) = T (A = j^t), 
which we have characterized in the previous section. 

The dynamical (four-point) correlation function is de- 
fined as the spatial correlation of the q field: 

G 4 (R,a,r) = (Sqj (a,r) 5q] {^r))\ t ^_^ =R (9) 

where 5q\ (a, r) = q\ (a, r) — Q (a, r) and is plotted 
for three packing fractions on fig. ©-left. Its sum 
over all R's defines the so-called dynamical susceptibil- 
ity X4 (a, r). Through simple manipulations, it is easy to 
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FIG. 6: Left: Four-point correlator G%(R) — 

G 4 (i?;cr (r*) ,t*) as a function of R for = 0.8402 (red), 
0.8417 (black), 0.8426 (blue). Right: Dynamical correlation 
as a function of the packing fraction; (Inset: Rescaling of 
\og[Gt(R)/Gt(0)] as a function of yj R/£,4 for 8 densities 
around <f>j). From 

show that X4 ( a -> T ) is related to the variance of Q as (see 
e.g. 0] for a review): 



X4 (a,r) = 




We have shown in [1] that X4 ( a -> T ) nas ? f° r a given 0, an 
absolute maximum xl — X4( a *? r *)? which as expected 
sits on the line corresponding to Q (a,r) w 1/2, i.e. such 
that half of the particles have moved by more than a. The 
amplitude of this maximum, which can be interpreted as 
a number of dynamically correlated grains, grows as (j) 
approached </>j, indicating that the jump motion of the 
grains becomes more and more collective as one enters the 
superdiffusive Levy regime, as anticipated above. Note 
in particular that r* behaves in the same way as r s jj , the 
time at which the superdiffusion exponent v reaches its 
maximum. We furthermore found [l| that the four-point 
correlation G%(R) = G±(R, «*,t*) is a scaling function 
of R/^4 (see inset of Fig. El-right), where £4(0) is the 
dynamical correlation length such that oc £f . 

All these results were reported in |1[ and are re- 
called here for completeness. The behavior of xtW 
was furthermore shown in [2] to be well accounted for 
by the following upper bound, derived in 0, Q: X4 > 
(dQ/d(j)) 2 (</> 2 ) c ? where (</> 2 ) c is the variance of the local 
density fluctuations. Here, we want to present further 
speculations, first on the relation between the size of the 
'micro-jumps' and dynamic correlations, and then on the 
higher moments of the distribution of Q (a,r). 

By conservation of the number of particles, the local 
change of density 5(j) is related to the divergence of the 
displacement field by: 5(j)/ = V • r. If we invoke a 
kind a Reynolds dilatancy criterion whereby the local 
density must fall below some threshold for the system 
to move, the displacement field must be correlated over 
some length £4 such that V • f ~ r/^4 w c, where c is a 
small number, possibly dependent on 0, and of the order 
of 10 -3 — 10 -2 (i.e. the relative difference between <pj 



and (j) g ). This immediately leads to a relation between 
the typical size of the jumps after time r and the required 
scale of the cooperative motion: 

&(t)~ — ocr^, r<r* (11) 
c 

This very simple argument predicts that £4 should be 
~ 10 2 — 10 3 times larger than typical displacements, 
which is indeed the case (see Fig. [6]). Furthermore, using 
r* oc td one finds r(r*) ex r max , and therefore, using 
the truncated Levy flight model above, a power-law di- 
vergence of £4(0) as \<j) — (j)j\~^ . This divergence might 
change if c strongly depends on the distance \<j) — <pj\. 

Turning now to higher cumulants of the distribution 
of Q (a*,r*), one can show that they are related to the 
space-integral of higher order correlation functions of the 
dynamical activity. For example, the skewness of Q is 
1/N 2 times the space integral of the 6-point correlation 
function, defined as: 

G 6 (R,R',a*,r*) = (Stfttfitf^-ZwiU-iL^ ' 

(12) 

where 8q\* is a shorthand notation for 5qj (a*,r*). 

Similar relations hold for higher moments. The sim- 
plest scenario is that all these higher order correlation 
functions are governed by the same dynamical correla- 
tion length £4, extracted from G^R^a^r). This is what 
happens in the vicinity of standard phase transitions, for 
example. If this is the case, and provided that xl oc £|, 
one can show that the following scaling relations should 
hold: 

/ ^4 X4 /1 o\ 

Q>~Q>,cY^y> ^8~^c,8-^ (13) 

where ^ , ^8 are respectively the skewness and kurtosis of 
1/N J2i Qi (related to the 6- and 8-point connected corre- 
lation functions), and c; Cj 6,^c,8 the corresponding values 
at the critical point, i.e. for systems of size TV smaller 
than the correlation volume xl oc £f. If these scaling 
relations are valid, the determination of c; C) 6,ft c ,8 using 
the above equations should give the same values for any 
(j) close to (j)j. Unfortunately, our statistics is not suffi- 
cient to make definitive statements, although the data is 
indeed compatible with such scalings. The notable fea- 
ture is that both c; C) 6 — —1 and c± —5 are found to 
be negative, although the error bar on both quantities 
is large. Interestingly, if we assume that the displace- 
ments are perfectly correlated within a correlation blob 
of size £4, the Levy flight model with \i = 1.6 makes the 
following predictions: 

^ c ,6 w -0.12; k; c , 8 w -1.37, (14) 

in qualitative agreement with our data (note however 
that there is an unknown proportionality factor in 
eq. (Q3J). A Gaussian diffusion model, on the other hand, 
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predicts c Cj e = and = —6/5, corresponding to a 
uniform distribution of Q in [0,1]. A kurtosis smaller 
than —6/5 means that the distribution of Q tends to be 
spiked around and 1 . Finally, as (j) increases beyond (j) j , 
the negative skewness of Q markedly increases. This is 
a sign that the dynamics becomes more and more inter- 
mittent, "activated", with a few rare events decorrelat- 
ing the system completely, while most events decorrelate 
only weakly. Since we fix Q (a*, r*) = 1/2, this leads to a 
diverging negative skewness in the limit where the proba- 
bility of rare decorrelating events tends to zero. The idea 
of characterizing the skewness of Q might actually be an 
interesting tool to characterize the strength of activated 
events in other glassy systems. 

CONCLUSION 

The central finding of this work is that the superdiffu- 
sive motion of our frictional grains in the vicinity of the 
jamming transition appears to be a genuine Levy flight, 
but with 'jumps' taking place on a Lilliputian scale. 
This interpretation contrasts with the original suggestion 
made in |l(, where anomalous diffusion was ascribed to 
persistent, spatially correlated currents. The vibration 
of the plate induces a broad distribution of jumps that 
are random in time, but correlated in space, and that 
can be interpreted as micro-crack events on all scales. 
As the volume fraction departs from the critical jam- 
ming density </>j, this distribution of jumps is truncated 
at smaller and smaller jump sizes, inducing a crossover 
towards standard diffusive motion at long times. This 
picture severely undermines the usefulness of harmonic 
modes as a way to rationalize the dynamics of our sys- 
tem (although this conclusion might of course not carry 
over to frictionless grains, or to thermal systems). The 
detailed study of these modes, and the difficulty to ana- 
lyze them in the present system, is discussed in [l5| . 

We have also presented several other speculations 
about the relation between the dynamical correlation 
length £4 and the size of the jumps, and the structure 
of higher order dynamical cumulants. The idea of using 
the 6-point skewness as a quantitative measure of the 
importance of activated events in the dynamics of glassy 
systems seems to us worth pursuing further. 

Note added and acknowledgment: We want to thank 
L. Berthier for interesting discussions about a preprint 
of his and collaborators [16| (that appeared in the very 
last stages of the present work), and where effects simi- 
lar to those discussed here are observed numerically in a 



model system. We thank him in particular for insisting 
on the importance of measuring velocity correlations in 
our system. 
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